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Abstract 

In about 95% of cases, congenital adrenal hyperplasia (CAH) is caused by mutations in 
CYP21A2 gene encoding steroid 21-hydroxylase (210H). Recently, we have reported four 
novel CYP21A2 variants in the Norwegian population of patients with CAH, of which 
p.L388R and p.E140K were associated with salt wasting (SW), p.P45L with simple virilising 
(SV) and p.V21 1 M + p.V281 L with SV to non-classical (NC) phenotypes. We aimed to 
characterise the novel variants functionally utilising a newly designed in vitro assay of 210H 
enzyme activity and structural simulations and compare the results with clinical phenotypes. 
CYP21A2 mutations and variants were expressed in vitro. Enzyme activity was assayed by 
assessing the conversion of 17-hydroxyprogesterone to 1 1-deoxycortisol by liquid 
chromatography tandem mass spectroscopy. PyMOL 1 .3 was used for structural simulations, 
and PolyPhen2 and PROVEAN for predicting the severity of the mutants. The CYP21A2 
mutants, p.L388R and p.E140K, exhibited 1.1 and 11.3% of wt 210H enzyme activity, 
respectively, in vitro. We could not detect any functional deficiency of the p.P45L variant 
in vitro; although prediction tools suggest p.P45L to be pathogenic. p.V21 1M displayed 
enzyme activity equivalent to the wt in vitro, which was supported by in silico analyses. 
We found good correlations between phenotype and the in vitro enzyme activities of the SW 
mutants, but not for the SV p.P45L variant. p.V21 1 M might have a synergistic effect together 
with p.V281L, explaining a phenotype between SVand NC CAH. 
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Introduction 

The CYP21A2 gene encodes the enzyme steroid 
21-hydroxylase (210H), which is essential for steroid 
synthesis in the adrenal cortex. Mutations in CYP21A2 
are the main cause of the autosomal recessive disorder 
congenital adrenal hyperplasia (CAH) (1). The classical 
forms of CAH are the salt wasting (SW) and simple 



virilising (SV) phenotypes with adrenocorticotropic 
hormone-driven excess of androgens leading to ambi- 
guous genitalia in female newborns. The SW phenotype 
has complete or almost complete loss of enzyme activity, 
leading to lack of Cortisol and aldosterone, whereas the 
SV displays various degrees of reduced enzyme activity 
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resulting in Cortisol insufficiency (1, 2). The non-classical 
(NC) CAH is milder and characterised by signs of 
hyperandrogenism postnatally and in adulthood, and 
associated with minor mutations (3). 

The most common mutations in the CYP21A2 gene are 
derived from a non-functional pseudogene CYP21A1P (4). 
Both CYP21A2 and CYP21A1P are located in the RCCX 
module of the HLA class III locus of chromosome 6 and share 
about 98% sequence homology. In general, there is a high 
correlation between the CYP21A2 genotype and the CAH 
phenotype (5, 6). In addition to the pseudogene mutations, 
more than 160 other pathological variants of CYP21A2 
have been reported (www.cypalleles.ki.se/cyp21.htm). 

In the present study, we used a new in vitro assay of 
21 OH enzyme activity and molecular computer modelling 
to predict the phenotypes of recently detected variants 
(p.L388R p.E140K, p.P45L and p.V211M) in the Norwegian 
population (7). We also checked a novel variant (p.A159T), 
detected in Norwegian patients with autoimmune 
Addison's disease (AAD) and in control populations. 

Subjects and methods 

Subjects 

We detected four novel variants: p.L388R, p.E140K, 
p.P45L and p.V211M, in a cohort of 59 patients diagnosed 
with CAH (7). Four unrelated female patients, each had 
one of the missense variants listed above, except the 
subject with p.V211M who also had a coupled p.V281L 
mutation on the same allele. All patients exhibited 
a deletion in the other allele (7). Clinical characteristics 
of the patients are listed in Table 1. The fifth variant, 
p.A159T, was found in heterozygosity in three of 370 
patients with AAD and in one out of 330 Norwegian 
healthy subjects. 

Genetic analyses of CYP21A2 

Protocols for sequencing and copy number analyses of the 
CYP21A2 gene were described previously (7, 8, 9, 10, 11). 

Mutagenesis 

Stratagene's QuikChange II Site Directed Mutagenesis Kit 
(Agilent Technologies, Santa Clara, CA, USA) was used to 
introduce point mutations in the pcDNA2 plasmid 
containing CYP21A2 cDNA (a kind gift from Prof. Holger 
Luthman, Department of Molecular Medicine, Karolinska 
Institute, Stockholm, Sweden (12)), corresponding to the 



Table 1 Clinical description of the CAH patients with novel 
variants in CYP21A2. 







Age of 




Variant 


Clinical phenotype 


diagnoses 


Treatment 


p.L388R 


Salt wasting 
Clitoris hypertrophy 


Newborn 


Mineralo- 

corticoids 
Glucocorticoids 


p.E140K 


Salt wasting 


1 week 
old 


Mineralo- 

corticoids 
Glucocorticoids 


p.P45L 


Simple virilising 
Clitoris hypertrophy 


4 years 


Glucocorticoids 


p.V211M + 


Simple virilising/ 


21 years 


Glucocorticoids 


p.V281L 


non-classical CAH 
Pubes growth and 
pronounced 
hirsutism from 
7 years old 







novel CYP21A2 variants (p.L388R, p.E140K, p.P45L, 
p.V211M and p.A159T). Three control variants of 
CYP21A2 were also generated: p.R102K (normal variant), 
p.I172N (SV phenotype) and p.G291R (SW phenotype). 
Primers used for mutagenesis are listed in Table 2. The 
mutagenesis was carried out according to the manufac- 
tures protocol, using 16 cycles with 5 min extension time. 
The mutagenesis was verified with direct sequencing, and 
plasmids from positive clones were purified by HiSpeed 
Plasmid Maxi Kit (Qiagen). 

In vitro protein expression of 21 OH mutants 

Non-radioactive TNT Quick Coupled Transcription/Tran- 
slation system (Promega) with Sp6 RNA polymerase was 
used for in vitro protein expression. Expression of all the 
210H variants was verified by western blot using a primary 
antibody against 210H (cl7, goat polyclonal IgG, 
Santa Cruz Biotechnology) and alkaline phosphatase- 
conjugated donkey anti-goat IgG secondary antibody 
(Santa Cruz Biotechnology). 

Relative quantitation of 21 OH from in vitro expression 

ELISA plates (MaxiSorp, Nunc Thermo Scientific, Wal- 
tham, MA, USA) were coated over night with 10 jj.1 of the 
in vitro expression product diluted in PBS, pH 7.4. Primary 
antibody against 210H (Santa Cruz Biotechnology), 
together with the corresponding alkaline phosphatase- 
conjugated secondary antibody (Sigma) and the FAST 
pNPP detection system (Sigma), was used to determine the 
relative amount of 2!OH in the reaction products. 
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Table 2 Primers used in mutagenesis. 



Mutations Sequences (5 — 3 ) 



p.P45L F: CAGCCCGACCTCCTGATCTATCTGCTTGG 

R: CCAAG C AG ATAG ATCAG GAG GTCG G GCTG 

p.R102K F: CAAGCTGGTGTCTAAGAACTACCCGGAC 
R: GTCCGGGTAGTTCTTAGACACCAGCTTG 

p.E140K F: CATGGAGCCAGTGGTGAAACAGCTGACCCAG- 
GAG 

R: CTCCTGG GTCAG CTGTTTCACCACTG G CTCCATG 
p.A159T F: GGCACCCCTGTGACCATTGAGGAGG 

R: CCTCCTCAATG GTCACAGG G GTGCC 
P.I172N F: CACCTGCAGCATCAACTGTTACCTCACC 

R: GGTG AG GTAACAGTTG ATGCTG CAG GTG 
p.G291R F: CTGCAGTGGACCTCCTGATCAGGGGCACTGAGAC- 
CACAGCAAAC 

R: GTTTGCTGTGGTCTCAGTGCCCCTGATCAG- 
GAGGTCCACTGCAG 
p.L388R F: GTCATCATTCCGAACAGACAAGGCGCCCACCTG 

R: CAG GTG G G CG CCTTGTCTGTTCG G AATG ATG AC 



The ELISA method was verified by selected reaction 
monitoring (SRM) against two signature peptides of 21 OH 
(p.85-98, LQEELDHELGPGASSSR and p.3 17-333, WAD- 
FAGRPEPLTYK) . Briefly, 5 nl of each version of 21 OH 
reaction products were taken out for quantitative 
determination of 210H. The lysates were purified on 
SDS-PAGE gel, and the bands corresponding to the 
molecular weight of wt 21 OH (56 kDa) were cut out for 
trypsin treatment. Stable isotope-labelled synthetic 
peptides (N-terminal replacement of lysine ( 13 C6, 15 N 2 ) 
or arginine ( 13 C 6 , 15 N 4 ) respectively) were added in 
constant amount before purification, resulting in a mass 
difference of 8 or 10 Da to the corresponding endogenous 
peptide. A Q-Trap 5500 (AB SCIEX) connected to a Dionex 
Ultimate NCR-3500RS LC system was used for the SRM 
analyses. The SRM data were analysed using Skyline, 
version 2.1 (https://skyline.gs.washington.edu/), and the 
most abundant transition was used for quantification. The 
SRM analyses were carried out using standard procedures 
by the Proteomics Unit (PROBE) at the University of 
Bergen. For quantitation, the amounts of 210H constructs 
were calculated relative to the amounts of the 21 OH 
wt protein. 

Enzyme activity assay 

The activity of the different 210H constructs was estimated 
by quantitating 11-deoxycortisol (11DOC) using 
17-hydroxyprogesterone (170HP) as substrate. The 
enzyme reactions were set up mainly as described 
previously (13). Briefly, from each 210H construct, 50 |il 
of the in vitro expression product were mixed with 15 pmol 



of recombinant human cytochrome P450 oxidoreductase 
(CYPOR, Sigma-Aldrich) and 25 ul of 300 uM 170HP (in 
10% methanol) on ice. Pre-warmed reaction buffer (50 mM 
potassium phosphate (pH 7.4), 20% (v/v) glycerol, 2 mM 
dithiothreitol, 0.1 mM EDTA and 0.0034% (v/v) Emulgen 
913 (Karlan Research Products Corporation, Cottonwood, 
AZ, USA)) was then added to a total volume of 500 liI. After 
60 s of pre-incubation at 37 °C, the enzymatic reaction was 
initiated by adding 5 |.d of 100 mM NAPDH. The enzymatic 
reaction was carried out for lh at 37 °C, and then 
terminated by putting onto ice. Further, the enzyme 
reaction solution was diluted 1:10 or 1:100 in a solution 
of 50% methanol and 0.1% formic acid in double-distilled 
(dd) H 2 0 (v/v). An internal standard (IS) solution (50 nM 
170HP-d8 and 50 nM HDOC-d2 in methanol) was then 
added to the diluted samples in a 1:1 ratio. Working 
standards of 11DOC in the range of 0.61-150 nM were 
prepared by diluting stock solutions in a solution of 50% 
methanol and 0.1% formic acid in ddH 2 0 (v/v) and adding 
IS (1:1). Working standards and 11DOC produced from 
each enzymatic reaction were quantified by liquid 
chromatography coupled to tandem mass spectroscopy 
(LC-MSMS) (14). The enzyme activity was calculated as the 
% ratio of 11DOC production between the 210H mutants 
to the wt, adjusted for the quantity of 210H wt and mutant 
in the enzymatic reactions. 

Structural simulations and in silico analyses 

The molecular graphic images were generated from a 
humanised version of bovine 210H as template (15), 
kindly provided from Dr Shozeb Haider (Centre of Cancer 
Research and Cell Biology, Queens's University of Belfast, 
Belfast, UK), using The PyMOL Molecular Graphics 
System, version 1.3 (Schrodinger, LLC, Portland, OR, 
USA). Variants were modelled using the Mutagenesis tool 
in PyMOL and distances between molecules measured 
using the Measurement tool. The severity of the variants 
was predicted by PolyPhen-2 v2.2.2r398 HumVar model 
(http://genetics.bwh.harvard.edu/pph2/) (16) and PRO- 
VEAN (http://provean.jcvi.org/index.php) (17) using 
default settings. Sequence alignments of 210H from 
seven different species (GenBank: AAB59440.1 {Homo 
sapiens), AAA30487.1 (Bos taurus), ABC69211.1 (Felis 
catus), AAD05573.1 (Rattus norvegicus) and BAF63009.1 
(Gallus gallus); NCBI reference sequence: XP_003311237.1 
(Pan troglodytes) and NP_999598.1 (Sus scrofa)) were carried 
out by default settings using Clustal omega (http://www. 
ebi.ac.uk/Tools/msa/clustalo/) (18) in order to assess the 
evolutionary conservation of the mutated amino acids. 
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Statistical analyses 

Statistical analyses were performed by the GraphPad Prism, 
version 5.02 Software (La Jolla, CA, USA). Correlations 
between the ELISA and SRM methods for quantifying the 
210H proteins produced were calculated by Pearson's 
correlation test. Deviation of enzyme activity of the 210H 
mutants compared with the wt 21 OH was calculated by 
ANOVA combined with Dunnett's multiple comparison 
tests using a significance level of P< 0.05. 

Results 

Verification and relative quantitation of 210H by 
ELISA and SRM methods 

Western blot analysis of the samples from in vitro expression 
showed a 56-kDa band corresponding to the molecular 
weight of 21 OH for all constructs (data not shown). 

The ELISA method for quantitation of the mutated 
21 OH proteins correlated very well with the SRM method 
(P<0.0001, Pearson's r=0.9206, and 95% CI (0.8016, 
0.9695), h = 19). Thus, the ELISA method was used for 
further quantitation of 21 OH. 

In vitro 21 OH activity studies 

The assay showed no significant reduction in the enzyme 
activity of the p.R102K normal variant relative to the wt 
(119.7%), whereas the p.I172N (SV phenotype) and 
p.G291R (SW phenotype) mutations showed significant 
(P<0.05) decline (1.6 and 0%, respectively, relative to wt 
(Table 3)). The novel CAH mutations, p.L388R and 
p.E140K, displayed significantly (P<0.05) reduced activity 



of 1.1 and 11.3% compared with wt respectively (Table 3). 
The other two novel CAH variants, p.P45L and p.V211M, 
were similar to the wt (105 and 99.5% respectively) as was 
p.A159T (126.6%) (Table 3). 

Structural simulations 

The structures of 21 OH and the novel variants based on the 
humanised model (15) are illustrated in Fig. 1. For p.L388R, 
the shift from the hydrophobic leucine to a negatively 
charged arginine in position 388 will probably lead to a 
disruption in a hydrophobic pocket surrounding the heme 
group (Fig. la). This is a highly conserved region (Fig. 2), and 
the p.L388R mutation was predicted as pathologic (Table 3). 

The p.E140 residue forms a salt bridge with the R444 
residue, and the shift from negatively charged glutamic 
acid to positively charged lysine breaks the salt bridge 
(Fig. lb). A negatively charged amino acid in this residue is 
highly conserved (Fig. 2); however, PolyPhen2 predicted 
this mutation as possible damaging, while PROVEAN 
predicted it as neutral (Table 3) in line with significant 
residual enzyme activity (see above). 

p.P45 is located early in the A helix (Fig. lc), which is a 
highly conserved region (Fig. 2). The proline-to-leucine 
shift in position 45 was predicted as pathologic by both 
PolyPhen2 and PROVEAN (Table 3), but in vitro enzyme 
activity was normal. 

The p.V211 residue (Fig. Id) may interact with L65 
through hydrophobic interactions and is located in a 
semiconservative region (Fig. 2). The shift from valine to 
methionine was predicted as benign/neutral (Table 3) in 
line with normal enzyme activity. 

The variant p.A159T (Fig. le) is located on the surface 
of the 210H protein. Sequence alignment shows that the 



Table 3 Enzyme activity (% of WT) from four to five different experiments and in silico prediction of CYP21A2 mutations and 
variants. 



Mutant/variant 


% Average enzyme 
activity (s.d.) 


Phenotype 


PolyPhen2 


Predicted phenotype (score) 
PROVEAN" 


p.K102R 


119.7 (22.5) 


Normal 


Benign (0.001) 




Neutral (0.722) 


p.A159T 


126.6 (29.9) 


Unknown 


Benign (0.023) 




Neutral (0.36) 


p.V211M 


99.5 (32.4) 


SV/NC b 


Benign (0.029) 




Neutral (h-1.019) 


P.I172N 


1.6 (0.8) 


SV 


Probably damag 


ng (0.998) 


Deleterious (h- 5.499) 


p.P45L 


105.0 (10.6) 


SV 


Probably damag 


ng (1.0) 


Deleterious (h- 5.067) 


p.G291R 


0 C 


SW 


Probably damag 


ng (1.0) 


Deleterious (h-4.637) 


p.E140K 


11.3 (2.4) 


SW 


Possibly damaging (0.49) 


Neutral (h- 1.208) 


p.L388R 


1.1 (0.6) 


SW 


Probably damag 


ng (0.986) 


Deleterious (h-4.995) 



NC, non-classical CAH; SV, simple virilising; SW, salt wasting. 

"Threshold as deleterious set to < 2.5. 

b Phenotype of p.V21 1 M + p.V281 L/deletion genotype. 

"Production of 11 DOC concentration under the detection limit of LC-MSMS method. 
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Figure 1 

Tertiary structure of 21-hydoxylase with positions of five novel variants 
generated by PyMOL v.3.1 . (a) p.L388R, (b) p.E140K, (c) p.P45L, (d) p.V211M 
and (e) p.A159T. 



alanine-to-threonine shift in position 159 is equal to the 
sequences of B. taurus and 5. scrofa 210H (Fig. 2), and was 
predicted as harmless (Table 3), again in line with a normal 
enzyme activity. 

Discussion 

In the present study, we introduced a novel method for 
functional studies of CYP21A2 mutants. Using the 
coupled transcription/translation system for protein 
expression, we were able to measure the 210H enzyme 
activity without the need of cell culture facilities and 
radiolabelling substrates, which is the current standard 
method (19, 20, 21, 22). We then characterised five novel 
and three known variants of the CYP21A2 gene found in 
patients with CAH (7) and AAD. Overall, the measured 
enzyme activities and structural simulation corresponded 
to the phenotype with some exceptions. 

The three polymorphisms used for verification of the 
in vitro enzyme assay behaved as expected; p.R102K (23) 
showed no reduction in activity, whereas the well-known 
SV mutation p.I172N had <2% activity and the G291R 
showed no activity at all, in accordance with previously 
reports (24, 25, 26). The novel p.L388R mutation showed 
almost complete abolishment of enzyme activity in vitro, 
which is in total agreement with the SW phenotype and 
the structural simulation. p.L388 is a crucial residue as it 
interacts with the heme prosthetic group through hydro- 
phobic interactions. When the leucine at this position 
is replaced by the positively charged and more space- 
demanding arginine, the heme group is restricted from its 
normal binding to 21 OH, which would have deleterious 
consequences for the enzymatic activity as well as the 
three-dimensional folding of 210H. 

The patient with the p.E140K mutation was diagnosed 
as SW, but even if in vitro enzyme assay showed a 
significant reduction in 210H activity (11%), it was still 
much higher than expected for the SW phenotype. In 
addition, the prediction tools did not foresee a severe 
damaging effect, and was more in line with a NC 
phenotype. However, structural analysis shows that the 
p.E140K disrupts a salt bridge with p.R444, which is crucial 
for maintaining the tertiary structure (15) that could 
explain the SW phenotype. Thus, prediction tools alone 
are not always precis (27). On the other hand, phenotype 
and genotype divergence in classical CAH has been 
described, and the same mutation can give different 
phenotypes even in siblings (28). 

The variant p.P45L was detected in a SV patient, but the 
in vitro enzyme activity did not correlate to such a severe 
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p.P45 



p.E140 



p.A159 



p.V211 



P.L388R 



Homo sapiens 
Pan troglodytes 
Bos taurus 
Sus scrofa 
Felis catus 
Rattus norvegicus 
Gallus gallus 



LQPDLPI YLLG ... 

LQPDLPI YLLG ... 

LQPNLPIHLLS ... 

LQPNLPI YLLG ... 

LQPDLPI YLLG ... 

LQPNLPVYLFG ... 

LHPGGPRHLHI ... 



ME PWEQLTQE 
MEPWEQLTQE 
MEPWVDQLTQE 
MEPRVEQLTQE 
MEPLVEQLTQE 
MEPLVEQLTQE 
LEPLLELQGWA 



PGTPVAIEEEF 
PSTPVAIEEEF 
AGAPVTIQKEF 
AGTPVTIQKEF 
AGTPVAIQKEF 
AGASVAIHKEF 
GAAPVDLFEAF 



IQIVDVIPFLR ... 

IQIVDVIPFLR ... 

IQILDMVPFLR ... 

IQILDMVPFLR ... 

IQVLDIVPFLR ... 

VQILDI I PFLR ... 

VRALDVLPSLR ... 



VI I PNLQGAHL 
VI I PNLQGAHL 
WI PNLQGAHL 
WI PNLQGAHL 
VVI PNLQGAHL 
I I IPNIQGANL 
I VI PNLFAAHH 



Figure 2 

Amino acid sequence alignments of 21-hydroxylase for the p.P45, p.E140, p.A1 59, p.V21 1 and p.L388 residues for seven species. 



phenotype. However, this variant was predicted as probably 
damaging/deleterious by computer modelling. The p.P45L 
residue is located in the N-terminal region of the enzyme 
close to the hydrophobic domain that anchors 210H to the 
endoplasmic reticulum (ER) membrane (15). The import- 
ance of the far N-terminal region has been experimentally 
determined earlier with regards to membrane targeting and 
insertion, as well as in vivo protein stability (29). It has also 
been demonstrated that the first 20 amino acids are 
necessary for membrane integration but that enzymatic 
activity in vitro can still be retained when these residues are 
deleted (30). Moreover, the p.P45 residue marks the start of 
an alpha helix that bends the direction of the polypeptide 
chain back towards the membrane. It is possible that this 
bending, for which proline in the 45 position is mandatory, 
is crucial for the correct orientation of 21 OH in the ER 
membrane. If p.P45 is replaced by another amino acid, the 
overall stability and in vivo enzymatic activity of the protein 
may be severely decreased. Another possibility is that the 
removal of the proline residue at position 45 alters the 
enzymes orientation to the ER membrane and somehow 
disturbs the electron transfer mediated in vivo by the 
membrane-bound CYPOR. As we added exogenous CYPOR 
in our in vitro assay, this effect could be lost. However, it is 
still possible that the rest of the enzyme maintains its 
structural integrity and that the intrinsic enzymatic 
function is intact. Similar effects of mutated proline residues 
have been described for other proteins, such as cubilin (the 
intrinsic factor-BT2 receptor) or the guanylate cyclase- 
activating protein 1 (involved in central vision) (31, 32). 
As we have carried out the functional characterisation of the 
p.P45L mutant in the absence of membranes, the severely 
reduced in vivo activity relative to wt 21 OH could be missed 
in our system. We therefore propose that future studies 
should compare our novel method using cell-free expression 
of 21 OH with the traditional use of transiently transfected 
COS 7 cells, preferably with LC-MSMS-based quantification 
of steroids. In addition to the p.P45L mutant, such a 



comparative study should include a wide range of known 
and well-characterised 210H mutations and variants. 

The patient with the p.V211M variant was classified 
with a phenotype between SV and NC CAH. She also carried 
the NC CAH mutation p.V281L on the same allele, whereas 
the second allele was deleted. The in vitro 210H activity of 
p.V211M was similar to the wt; thus the cause of the 
disorder might be from the p.V281L mutation, which 
actually has been reported with SV phenotype (6). An 
almost similar mutation, p.V211L, has been described 
previously in a patient with NC CAH together with 
p.V281L, and was suggested as normal based on no 
conservation between four mammalian species (33). 
Structural analysis shows that p.V211 is a part of the 
substrate cavity site and that it probably interacts with L65 
through hydrophobic interactions (15). The conversion of 
valine into methionine at this position may alter this 
interaction slightly, but probably not to a damaging degree. 
PolyPhen2 and PROVEAN predicted p.V211M as benign/ 
neutral, but with a score close to the threshold. It is unclear 
whether this mutation contributes to the intermediate form 
of SV/NC CAH due to a synergistic effect, similarly that has 
been seen in other studies were mild NC variants together 
have given SV (20, 21, 22, 34, 35) or SW phenotypes (36). 

Finally, the p.A159T variant showed no reduction in 
the enzyme activity in vitro. This variant was found in 
heterozygosity in three patients with AAD and in one 
healthy subject, and was predicted as harmless. The 
p.A159 residue is located on the surface of the protein in 
a region outside the secondary structure motifs and away 
from the active site, and occurs naturally in other 
mammals such as pigs and cows (Fig. 2). The substitution 
of alanine with threonine is also common in homologous 
proteins from related species. 

In conclusion, the in vitro enzyme assay and structural 
simulations revealed good correlations to the most severe 
phenotypes. For the other variants, a heterogeneous 
picture was observed. Nevertheless, in vitro enzyme assays 
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and structural simulations are valuable supplements for 
evaluating the severity of novel CYP21A2 variants. 
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